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Исследованы собственные движения пилы путём построения и анализа математической модели инст! румента 
на основе уравнений теории упругости. Предложенная авторами математическая модель может быть ис- 
пользована при построении систем диагностики протекания технологического процесса. 
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Введение. Повышенная вибрация и шумоактивность при работе пилорамы вызывает необходи- 
мость в особых условиях закрепления полотна, что позволило бы изменить амплитудно-частотные 
характеристики пилы с целью отстройки собственных частот от частот возбуждений, определяе- 
мых рабочим процессом. С этой целью проведём исследование собственных движений пилы при 
следующих условиях. 

Результаты исследований. Считаем, что две боковые стороны пластины упруго защемлены, а 
две другие свободны, уравнение колебаний запишем в виде: 
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где р — толщина пластины; р — плотность стали; р Ри) — цилиндрическая жесткость 
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пластины; //, — сила предварительного натяжения пилы [1]. 


Продольные силы предварительного натяжения изменяют частоты её собственных коле- 
баний. Кроме того особенности граничных условий определяются также наличием «упругого за- 
щемления». 

Граничные условия: 
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Таблица 1 
Характеристики пилы как объекта исследования 
Параметры рамной пилы Характеристики 
В— ширина пилы 160-180 
В — толщина пилы 1.6; 1.8; 2.0; 2.2; 2.5 
[ — длина пилы 1100; 125; 1600; 1750; 1950 
р — цилиндрическая жёсткость 75-275 
о — напряжение в серединной плоскости пилы | 80-120 МПа 
№, — сила натяжения пилы 20.5-450 кН 
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Эти данные будут использованы при формулировании необходимых допущений. 
В первом приближении, с учётом того, что соотношение ширины и длины пластины при- 


мерно В ‚ воспользуемся соотношением изменённой вследствие натяжения собственной часто- 


ты, приведённым для балки: 


Кг 
где «, — собственная частота для пластины при № =0; № — усилие предварительного натяже- 

ния; М, — критическая сжимающая сила для балки, соответствующая продольному изгибу в 

2 
п-Е7 
плоскости колебания /,, = в [2]. 
И Для указанных граничных условий можно принять решение 

уравнения в виде (2): 


и! (&,п,Ё) = С, с0$ К. . 0$ К.п -ехр (№) = Ат (&,п)-ехр (№), 


где ] — мнимая единица; К,,К, — волновые числа. 
а 
Вблизи границ & = о решение представляется в виде (3): 


и (Е, п,Ё) = оз К,п. [С, - с0$ К.Е + # (&) | -ехр (дё), 
(3) 


где Ё (=) — функция, позволяющая выполнить условия закрепле- 


ния. 
Ь 
Вблизи границ п = +5 решение представляется в виде: 


и (Е, п,Ё) = с0$ К. - [С - 60$ п + ф (п) | -ехр (№2), 
(4) 


где ф(п) — функция, позволяющая выполнить условия закрепле- 











НИЯ. 


Рис. 1. Расчётная схема пилы 


Следовательно, общее выражение Ат (п) имеет вид: 
Ат (&, п) = Со 0$ К.Е - 0$ К, + # (&)- с0$ КП + 4 (п) - с0$ АН . (5) 
В средней части пластины /(&) и 4(п) пренебрежимо малы, поэтому выражение 
Ат (&, п) = Со 0$ К, - с0$ К.П должно удовлетворять уравнению: 
У?У?Ат (&,п)- “Ат (&п)=0, (6) 
р.^ 


где У? — оператор Лапласа; а“ = о: 


Подставляя выражение Ат (п) в уравнение (1) получим: 
К! с0$ К.Е -С, с0$ КП + 2С, (-Кр соз А.Е) - (-К3 с0$ К.П) + К; с0$ Юй - Су с0$ К.Е — 
—а* с05 КЕ. с05КП-С, =0 
КА +2 К + а" =0, то есть а? =К?+К». (7) 
Вблизи границ &= = несущественно влияние функции ф(п) ‚ поэтому уравнение (6) 
можно записать в виде: 


УРУ" [1 (Е) созкол | - а* -^ (5)-с05№п=0. (8) 
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После дифференцирования имеем: 


41 (=) 4"! (&) 








= Юп+2——>- (-КЗ со$Кьп) + / (Е) - К соз п — а* 0$ КП =0 (9) 
в Ан +! (=) (К; -а*)=0. 


Корни соответствующего характеристического уравнения получаем в виде: а =К? +22. 
Выражение для функций / (&) получим в виде: 


Г (5) =С, ехр [-2 Е + 2] — вблизи границы & = г ‚ где а, = Ме +2К2 
Г (5) = С, ехр [-2 [2 =) — вблизи границы & -5, где а, = Ме +22 


Аналогично рассуждая, можно определить функцию ф(п). 





Й 





Вблизи границ п несущественно влияние функции г (5), поэтому уравнение (9) 


можно записать в виде: 
У?У? [ф(п)- 0$ КЕ |-а* -ф(п)-с0$ А.Е =0. (10) 
После дифференцирования имеем: 
а Я 
+ (п) 05 АЕ + аи. (-К? со$ К.Е )+ (п)! с0$ КЕ — а* с0$ Е = 0 
2 а? 
2 И ке). (6-6) =0. 
Корни соответствующего характеристического уравнения получаем в виде: а =2^2 +К?. 
Выражение для функций ф(п) получим в виде: 


+(п) =С, е[-0, |п 2] — вблизи границы п = =, где а, =| к? +? 
ф(п) = С, е[-0, [2- ") — вблизи границы п = : ‚ где а, = Маке +к 


Рассмотрим граничные условия на концах п = +7: 
0*Ат (Е, 0*Ат (Е, дАт (Е, 
Ат (&,п)=0 А ый О ЕН 
Г] 06 Г] 
где с — жёсткость упругого защемления (4); х — коэффициент Пуассона для стали. 











; (11) 








и $ Ь 
Учтём, что вдоль всей стороны пластины п = — функция / (&) равна нулю во всех точ- 


ках за исключением угловых точек, поэтому можно записать: 


Ат (Ел) [© сок, 2 (0) -о05КЕ, (12) 


Ат(Е, п] в = | 05 К, 5 „с, .СО5КЕ=О. 
2 


Для записи второго условия определим необходимые производные: 


аАт (5, 

ар с ил .2,64 0: [П- я] (13) 
а*Ат (Е, 

р = 05 КЕ. | с0$ КП — @2С, ехр |5 [п + 
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Подставляя выражения производных в уравнение (13), получим при П= -5: 


со$ К.Е. [- С.К? с0$ К, 2+0С ‚©, со$ К, 2+С, (-К? со А.Е) = 


= 5.03 КЕ [с 5 К, 5 = сс, 
Окончательно получаем систему уравнений вида: 


С, со5 К, 2+С, =0, 


Ь Обе аыниЫ С 
С [-в 0$ А, 5 КУ + т, эт, =] +6, |< — УК? +а, 5 =0. 


Рассмотрим условия свободных краёв: 
бАт(Еп)  „бАт(Еп) _‹ бАт(Е,п) Ат (Е,п) 


= > а 
6? оп? = 0Е3 ( ) Е0п? й 


при &=+— 

ри 5 5 
х и а 

Учтём, что вдоль всей стороны пластины & = в функция ф(п) равна нулю во всех точ- 


ках за исключением угловых точек, поэтому можно записать: 
Ат (п) = [С с05 А.Е + / (&)|-с05Кп =0, 


т (5.1) = | с05 А, 5+ с, -С05АП=0. 


Для записи второго условия определим необходимые мы 
а*Ат (Е, 
о = с0$ А.П. |- с0$ К, & — а?С, ео[-о [5 2] 
а?Ат (&, 
о = 60$ К.Е. [в с0$ К›П- а2С, ео[-в.[1 и] 
аЗАт (&,п) 
СЕ? 
а*Ат (5, п) 
аН* 
Граничные условия на свободных сторонах, учитывая, что косинус является чётной функ- 
цией, записываем в виде: 


= с05 А.П: си 51п А, & — а;С, ехр |-(: +2 "| 


=-К2 с05 К.П. ск эп А! — а, С, ехр |= | + =] 


с05 К.П: [- —С.К? с0$ К, — 2+0 |- УК? со С, С0$ К, — вы 0, 


05 К>П. [-< эп А, 5 _ ос.) + (2-м), с0$ ки - Со тк, — ОА а,С ‚)- 0. 
Окончательно получаем систему уравнений вида: 
а а 
С: |-& со5 А 5 - УК? 0$ К, = + С, (@-\к2)=0, 


е [бяка + -у)6зтк 2+ (а, (2-у) к -а)=0, 


где а, = \/К? +22 ; а, = 22 +К2. 
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Рис. 2. Зависимость функции Е (К, К,) от волновых чисел 


Для одновременного выполнения приведённых уравнений необходимо, чтобы определи- 
тели этих уравнений равнялись нулю, то есть решения следующих уравнений: 
Ь 
с0$ К, 5 ы 
Е1 = В Ь Ь =0 , 
ИЕ [о 
-К? со5 К, — - К2\с05 К, —+— К, $тК, — а? - УК? +а,.— 
2 р 2 1 2 2 р) 2 2 2 2 1 2. Б 


-К? 0$ К, 9 — УК? со, 9 а — УК? 
Е? = - 2 -0. 


К эп, 5 +(2-у) ик; 5 а, (2-\)К‚ — а 


Для обоснования дальнейших рассуждений был использован пакет прикладных программ 
МаЧаЬ. В результате была построена поверхность, определяемая выражением: 
Е(К,,К,) = аб$ (ЕК, ,К,))+ аЬ5 (Е2(К,,К,)) (рис. 2). 

Как следует из анализа полученной поверхности, «овраг» существует. Для того чтобы 
уточнить начальные значения волновых чисел, строим контурный график (рис. 3). Как видно из 
него, значения волновых чисел лежат около значений К, = 5, К, =5 (графики построены относи- 
тельно индексов, значит, при расчётном шаге 0,1 индексу 50 соответствует волновое число 5). 

Расчёт, проведённый в среде МаМаЬ с использованием встроенной функции «{о№е», по- 
зволил получить следующие результаты: 





К = Ка! = ех = 
5.194 8.3576е-015 1 
5.0265 -2.2737е-012 

















Заключение. Как видно из полученных результатов, значение ошибки решения задачи «Киа!» 
имеет порядок 10`'? -10-! , то есть точность вполне достаточна. Значение «ех» является положи- 
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тельным числом, что, по условиям МаЧаБ’а, означает успешное прохождение расчёта. Получен- 
ные результаты позволяют использовать разработанные программы для проведения математиче- 
ского моделирования как свободных, так и вынужденных движений пилы. 
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Рис. 3. Контурный график 
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